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Abstract 



Modifications to quark and antiquark fragmentation functions due to quark-quark 
(antiquark) double scattering in nuclear medium are studied systematically up to 
order 0(a 2 s ) in deeply inelastic scattering (DIS) off nuclear targets. At the order 
0{a1), twist-four contributions from quark-quark (antiquark) rescattering also ex- 
hibit the Landau-Pomeranchuck-Midgal (LPM) interference feature similar to gluon 
<^ \ bremsstrahlung induced by multiple parton scattering. Compared to quark-gluon 

scattering, the modification, which is dominated by t-channel quark-quark (anti- 
quark) scattering, is only smaller by a factor of Cf/Ca = 4/9 times the ratio of 
quark and gluon distributions in the medium. Such a modification is not negli- 
gible for realistic kinematics and finite medium size. The modifications to quark 
(antiquark) fragmentation functions from quark-antiquark annihilation processes 
are shown to be determined by the antiquark (quark) distribution density in the 
medium. The asymmetry in quark and antiquark distributions in nuclei will lead 
to different modifications of quark and antiquark fragmentation functions inside 
a nucleus, which qualitatively explains the experimentally observed flavor depen- 
dence of the leading hadron suppression in semi-inclusive DIS off nuclear targets. 
The quark-antiquark annihilation processes also mix quark and gluon fragmentation 
functions in the large fractional momentum region, leading to a flavor dependence 
of jet quenching in heavy-ion collisions. 
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1 Introduction 



Multiple parton scattering in a dense medium can be used as a useful tool to study 
properties of both hot and cold nuclear matter. The success of such an approach has been 
demonstrated by the discovery of strong jet quenching phenomena in central Au + Au 
collisions at the Relativistic Heavy- ion Collider (RHIC) [1,2,3] and their implications 
on the formation of a strongly coupled quark-gluon plasma at RHIC [4,5]. However, 
for a convincing phenomenological study of the existing and future experimental data, 
a unified description of all medium effects in hard processes involving nuclei, such as 
electron- nucleus (e + A), hadron- nucleus (h + A) and nucleus-nucleus collisions (A + 
A) has to be developed [6,7]. This must include the physics of transverse momentum 
broadening [8], strong nuclear enhancement in DIS [9] and Drell-Yan production [10,11], 
nuclear shadowing [12] , and parton energy loss due to gluon radiation induced by multiple 
scattering [13,14,15,16,17,18,19]. 

There exist many different frameworks in the literature to describe multiple scattering 
in a nuclear medium [20,21,22]. Among them the twist expansion approach is based on 
the generalized factorization in perturbative QCD as initially developed by Luo, Qiu 
and Sterman (LQS) [23]. In the LQS formalism, multiple scattering processes generally 
involve high-twist multiple-parton correlations in analogy to the parton distribution op- 
erators in leading twist processes. Though the corresponding higher twist corrections are 
suppressed by powers of 1/Q 2 , they are enhanced at least by a factor of A 1 / 3 due to 
multiple scattering in a large nucleus. This framework has been applied recently to study 
medium modification of the fragmentation functions as the leading parton propagates 
through the medium [18,19]. Because of the non-Abelian Landau- Pomeranchuck-Midgal 
interference in the gluon bremsstrahlung induced by multiple parton scattering in nuclei, 
the higher-twist nuclear modifications to the fragmentation functions are in fact enhanced 
by A 2 / 3 , quadratic in the nuclear size [18,19]. Phenomenological study of parton energy 
loss and nuclear modification of the fragmentation functions in cold nuclear matter [24] 
gives a good description of the nuclear modification of the leading hadron spectra in semi- 
inclusive deeply inelastic lepton-nucleus scattering observed by the HERMES experiment 
[25,26]. The same framework also gives a compelling explanation for the suppression of 
large transverse momentum hadrons discovered at RHIC [27]. 

The emphasis of recent studies of medium modification of fragmentation functions has 
been on radiative parton energy loss induced by multiple scattering with gluons. Such 
processes indeed are dominant relative to multiple scattering with quarks because of the 
abundance of soft gluons in either cold nuclei or hot dense matter produced in heavy-ion 
collisions. Since gluon bremsstrahlung induced by scattering with medium gluons is the 
same for quarks and anti-quarks, one also expects the energy loss and fragmentation 
modification to be identical for quarks and anti-quarks. However, in a medium with fi- 
nite baryon density such as cold nuclei and the forward region of heavy-ion collisions, the 
difference between quark and anti-quark distributions in the medium should lead to differ- 
ent energy loss and modified fragmentation functions for quarks and antiquarks through 
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quark- ant iquark annihilation processes. To study such an asymmetry, one must consider 
systematically all possible quark-quark and quark-antiquark scattering processes, which 
will be the focus of this paper. 

In this study we will calculate the modifications of quark and antiquark fragmentation 
functions (FF) due to quark-quark (antiquark) double scattering in a nuclear medium, 
working within the LQS framework for generalized factorization in perturbative QCD. 
For a complete description of nuclear modification of the single inclusive hadron spec- 
tra, one still have to consider medium modification of gluon fragmentation functions in 
addition to modified quark fragmentation function due to quark-gluon scattering [18]. 
The theoretical results presented in this paper will be a second step toward a complete 
description of medium modified fragmentation functions. However, one can already find 
that quark-quark (antiquark) double scattering will give different corrections to quark 
and antiquark FF, depending on antiquark and quark density of the medium, respec- 
tively. This difference between modified quark and antiquark FF may shed light on the 
interesting observation by the HERMES experiment [25,26] of a large difference between 
nuclear suppression of the leading proton and antiproton spectra in semi-inclusive DIS 
off large nuclei. Such a picture of quark-quark (antiquark) scattering can provide a com- 
peting mechanism for the experimentally observed phenomenon in addition to possible 
absorption of final state hadrons inside nuclear matter [28,29]. 

The paper is organized as follows. In the next section we will present the general for- 
malism of our calculation including the generalized factorization of twist-4 processes. 
In Section III we will illustrate the procedure of calculating the hard partonic parts of 
quark-quark double scattering in nuclei. In Section IV we will discuss the modifications 
to quark and antiquark fragmentation functions due to quark-quark (antiquark) dou- 
ble scattering in nuclei. In Section V, we will focus on the flavor dependent part of the 
medium modification to the quark FF's due to quark-antiquark annihilation and we will 
discuss the implications for the flavor dependence of the leading hadron spectra in both 
DIS off a nucleus and heavy-ion collisions. We will summarize our work in Section VI. In 
the Appendix A-l, we collect the complete results for the hard partonic parts for different 
cut diagrams of quark-quark (antiquark) double rescattering in nuclei. We also provide 
an alternative calculation of the hard parts of the central-cut diagrams in Appendix A-3 
through elastic quark-quark scattering or quark-antiquark annihilation as a cross check. 



2 General formalism 



In order to study quark and antiquark FF's in semi-inclusive deeply inelastic lepton- 
nucleus scattering, we consider the following processes, 

e{L l ) + A{p)-^e{L 2 ) + h{i h )+X , 
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Fig. 1. Lowest order and leading- twist contribution to semi-inclusive DIS. 



where Li and L 2 are the four momenta of the incoming and outgoing leptons, and lh 
is the observed hadron momentum. The differential cross section for the semi-inclusive 
process can be expressed as 



where p = [p + ,0,0±] is the momentum per nucleon in the nucleus, q = L 2 — L\ = 
[—Q 2 /2q~, q~, 0±] the momentum transfer carried by the virtual photon, s = (p + L\) 2 
the lepton-nucleon center-of-mass energy and «em is the electromagnetic (EM) coupling 
constant. The leptonic tensor is given by L MJ , = l/2Tr(7 • £i7 M 7 • L 2 'ju) while the semi- 
inclusive hadronic tensor is defined as, 



^ d*4 = 2 I>W0)|X, h)(X, h\J u (0)\A)2n5 4 (q +p -p x -£ h ) (2) 

where J2x runs over au possible final states and J M = J2 q £q4 ) qln 1 Pq ls the hadronic EM 
current. 

Assuming collinear factorization in the parton model, the leading-twist contribution to 
the semi-inclusive cross section can be factorized into a product of parton distributions, 
parton fragmentation functions and the partonic cross section. Including all leading log 
radiative corrections, the lowest order contribution [0(a£)] from a single hard 7* + q 
scattering, as illustrated in Fig. 1, can be written as 



dW, 



dz h 



/ dx f£{ x >tf) H $${ x >p><i) D <i->h(zh,tf l 



H$(x,p, q) = y Tr(7 ■ jrf^ ■ (q + xp)^ u ) -^~K X ~ x b) 



(3) 
(4) 
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where the momentum fraction carried by the hadron is defined as Zh = £^/q~, xb = 
Q 2 /2p + q~ is the Bjorken scaling variable, /i 2 - and fi 2 are the factorization scales for 
the initial quark distributions f q (x,n 2 ) in a nucleus and the fragmentation functions 
in vacuum D q _^h(zh, /i 2 ), respectively. The renormalized quark fragmentation function 
D q ^,h(zh, /x 2 ) satisfies the DGLAP QCD evolution equations [30]: 



dD q ^ h (z h , fi 2 



dD g ^ h (z h ,/j, 2 
d In fi 2 



a s (fi 2 ) f dz 

z 



2n 



lq^q g (z)Dg^h(z h /z,/J, 2 ) 

+ lq^gq{z)Dg^ h {z h /z,H 2 ) 



(5) 



a s (/i 2 ) f dz 

z 



2tt 



2nt 



Yl lg^qq{z)D q ^ h {z h / Z, /i 2 ) 
q=l 

+ -fg^gg(z)Dg_ th (Z h /Z,H 2 ) 



(6) 



where 'ja^bc(z) denotes the splitting functions of the corresponding radiative processes 
[31,32]. 

In DIS off a nuclear target, the propagating quark will experience additional scatterings 
with other partons from the nucleus. The rescatterings may induce additional parton 
(quark or gluon) radiation and cause the leading quark to lose energy. Such induced 
radiation will effectively give rise to additional terms in the DGLAP evolution equation 
leading to a modification of the fragmentation functions in a medium. These are power- 
suppressed higher-twist corrections and they involve higher-twist parton matrix elements. 
We will only consider those contributions that involve two-parton correlations from two 
different nucleons inside the nucleus. They are proportional to the thickness of the nucleus 
[18,23,33] and thus are enhanced by a nuclear factor A 1 / 3 as compared to two-parton 
correlations in a nucleon. As in previous studies [18,19], we will limit our study to such 
double scattering processes in a nuclear medium. These are twist-four processes and give 
leading contributions to the nuclear effects. The contributions of higher twist processes or 
contributions not enhanced by the nuclear medium will be neglected for the time being. 

When considering double scattering with nuclear enhancement, a very important process 
is quark-gluon double scattering as illustrated in Fig. 2. Such processes give the dominant 
contribution to the leading quark energy loss and have been studied in detailed in Refs. 
[18,19]. The modification to the vacuum quark fragmentation function from quark-gluon 
scattering is, 



+ 5{z 




(7) 
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Fig. 2. A typical diagram for quark-gluon double scattering with three possible cuts [central(C), 
left(L) and right(R)]. 

where the +function is defined as 

1 N 1 



dz- 



F(z) 



o 



dz 



o 



F(z)-F(l) 
l-z 



for any F(z) that is sufficiently smooth at z — 1 and the twist-four quark-gluon correla- 
tion function, 



Tf g (x,x L ) = J -^-dy^dy^e i{x+XL)p+y (1 - e ^ lXLP+ ^)(l - e~ lXLP ' 
x(A\j q (0) ^F a + (y^F + °(y^ q (y-)\A)6(-y^6(y- - ^ 



(y -Vi ) 



(9) 



has explicit interference included. The matrix element in the virtual correction [the term 
with 5(z — 1)] is defined as 



, 4) 



dz- 



1 - z 



2T£(x,x L )\ z=1 



[l + z 2 )T n A n ( Xj x L ) 



(10) 



Since T^(x,xl)/ f^(x) is proportional to gluon distribution and independent of the flavor 
of the leading quark, the suppression of the hadron spectrum caused by quark-gluon or 
antiquark-gluon scattering should be proportional to the gluon density of the medium 
and is identical for quark and antiquark fragmentation. It was shown in Ref. [24] that 
such modification of parton fragmentation functions by quark-gluon double scattering 
and gluon bremsstrahlung in a nuclear medium describes very well the recent HERMES 
data [25] on semi-inclusive DIS off nuclear targets. 
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Fig. 3. Diagram for leading order quark- antiquark annihilation with three possible cuts [cen- 
tral^), left(L) and right(R)]. 




Fig. 4. A typical diagram for next-to-leading order correction to quark-antiquark annihilation 
with three possible cuts [central(C), left(L) and right(R)]. 

In this paper, we will consider quark-quark (antiquark) double scattering such as the 
process shown in Fig. 3 and its radiative corrections at order 0{a 2 s ) in Fig. 4. The 
contributions of quark-quark double scattering is proportional to the quark density in 
a nucleon, while the contribution of quark-gluon double scattering is proportional to 
the gluon density in a nucleon; and the gluon density is generally larger than the quark 
density in a nucleon at small momentum fraction. However, as pointed out in earlier works 
[18], quark-quark double scattering mixes quark and gluon fragmentation functions and 
therefore gives rise to new nuclear effects. The annihilation processes as shown in Figs. 3 
and 4 will lead to different modifications of quark and antiquark fragmentation functions 
in a medium with finite baryon density (or valence quarks). Such differences will in 
turn lead to flavor dependence of the nuclear modification of leading hadron spectra as 
observed in HERMES experiment [25,26]. 

Quark-quark double scattering as well as quark-gluon double scattering are twist-4 pro- 
cesses. We will apply the same generalized factorization procedure for twist-4 processes 
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as developed by LQS [23] for semi-inclusive processes in DIS. In general, the twist-four 
contributions can be expressed as the convolution of partonic hard parts and two-parton 
correlation matrix elements. In this framework, contributions from double quark-quark 
scattering in any order of a s , e.g., the quark- ant iquark annihilation process as illustrated 
in Fig. 4, can be written in the following form, 



<WZ „ fp+dy 



dzh 



x(A\^^ q (y~)My^^ q (y^\A). (11) 

Here we have neglected transverse momenta of all quarks in the hard partonic part. 
Transverse momentum dependent contributions are higher twist and are suppressed by 
{k'jJ/Q 2 , Therefore, all quarks' momenta are assumed collinear, k% = x^p and k% = x^p. 

H av (y ~ i Hi 5 V2 i Vi <L z h ) is the Fourier transform of the partonic hard part H^ v (x,Xi,X2,p,q,Zh) 
in momentum space, 



H° u (y-,y^y^p,q,z h )= [ dx ^ ^ ^y-^p^+i^ ^) P ^ 



xH^(x,Xi,X2,p,q,z h ) 
dxH^J(x,p,q)H D (y~,y^,y^,x,p,q,z h ) , (12) 

where, in collinear approximation, the hard partonic part HjPj(x,p,q) [Eq. (4)] in the 
leading twist without multiple parton scattering can be factorized out of the high-twist 
hard part H® u (x, Xi, Xi,p, q, z^). The momentum fractions x, X\, and X2 are fixed by 
(^-functions of the on-shell conditions of the final state partons and poles of parton prop- 
agators in the partonic hard part. The phase factors in H ^ u (y~ , y±, y^iVi z h) can then 
be factored out, which in turn will be combined with the partonic fields in Eq. (11) 
to give twist-four partonic matrix elements or two-parton correlations. The quark-quark 
double scattering corrections in Eq. (11) can then be factorized as the convolution of 
fragmentation functions, twist-four partonic matrix elements and the partonic hard scat- 
tering cross sections. For scatterings (versus the annihilation) with quarks (antiquarks), 
a summation over the flavor of the secondary quarks (antiquarks) should be included 
in two-quark correlation matrix elements and both t, u channels and their interferences 
should be considered for scattering of identical quarks in the hard partonic parts. 

After factorization, we then define the twist-four correction to the leading twist quark 
fragmentation function in the same form [Eq. (3)], 



dW D 

r fj,i> 

dz h 



£ / dxf q A (x)H^(x,p,q)AD^ h (z h ) . (13) 
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3 Quark-quark double scattering processes 



In this section we will discuss the calculation of the hard part of quark-quark double 
scattering in detail. The lowest order process of quark-quark (antiquark) double scattering 
in nuclei is quark-antiquark annihilation (or quark-gluon conversion) as shown in Fig. 3. 
The hard partonic parts from the three cut diagrams in this figure are [18]: 



H o,c(V >fi '2/2 ,x,p,q,z h ) =D g ^ h ( Zh y^2C F ^e txp+r 



X0(-jfe)0(y--yr) 



2na s xb 



H o,L(y iyiiy2i x iV,qi z h)=D q ^ h {z h )—^-2C F —e ixp y 



x %i -y 2 ) d (y -Vi) 



H o,n(y iyi>y2,x,v,q,Zh) = D q ^ h {z h )—^-2C F —e lxp y 



27 ra s x B 

2Cf q~ 2 

x9(-V2) (V2 - Vi) ■ 



(14) 
(15) 
(16) 



The main focus of this paper is about contributions from the next-leading order correc- 
tions to the above lowest order process. There is a total of 12 diagrams for real corrections 
at one- loop level as illustrated in Fig. 5 to Fig. 16 in the Appendix A-l, each having up 
to three different cuts. In this section, we demonstrate the calculation of the hard parts 
from the quark-antiquark annihilation in Fig. 4 in detail as an example. We will list the 
complete results of all diagrams in Appendix A. 

One can write down the hard partonic part of the central-cut diagram of Fig. 4 (Fig. 5 
in Appendix A-l) according to the conventional Feynman rule, 



dz 



Zh< 



D g ^ h {—) / dx 



dx\ dx2 
2tt 27t' 



,iX!p-<-y +ix 2 p^y t 



i( x - Xl - X2 )p+y 2 r J^L T 

J (2tt) 4 



27T^ + (£ 2 )27T<5 + (^)(5(1 - z- 



jr _ C F 4 7 • (g + %lP) 7 • (g + XlP ~ I) i>_ 

(q + xip) 2 — ie ^ a (q + x\p — £) 2 — ie ^ 2 ^ a 



x 



7 ■ (q + xp — 



rU, , T ; (g t,T ) - ^(0^(O 



(q + xp — I) 2 + ie (q + xp) 2 + ie 



(17) 



where S + is a Dirac cfe/ta-function with only the positive solution in its functional variable, 
e ap (£) = —g ap + (n a £ p + n p £ a )/n ■ £ is the polarization tensor of a gluon propagator in 
an axial gauge (n ■ A = 0) with n = [1, 0~, 0j_], £ and £ g = q + {x\ + x 2 )p — £ are the 
4-momenta carried by the two final gluons respectively. The fragmenting gluon carries 
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a fraction, z = £ g jq , of the initial quark's longitudinal momentum (the large minus 
component). 

To simplify the calculation in the case of small transverse momentum ij- <C q~~,p + , we 
can apply the collinear approximation to complete the trace of the product of 7-matrices, 



H « 7 • 



1 

! 4t 



-Tr 



1-taH 



According to the convention in Eqs. (11) and (12), contributions from quark-quark double 
scattering in the nuclear medium to the semi-inclusive hadronic tensor in DIS off a nucleus 
can be expressed in the general factorized form: 



w qq,nv j2 J dxH$(x,p, q) J P ^ dy x dy 2 H°(y ,y 1 ,y 2 , x,p, q, z h ) 

x (A\^(0)^(y-)^ q (y^^)\A) . (19) 

After carrying out the momentum integration in x, Xi, X2 and £ ± in Eq. (17) with the 
help of contour integration and 5-functions, one obtains the hard partonic part, H° , of 
the rescattering for the central-cut diagram in Fig. 4 (Fig. 5) as 



dz h 



2(1 + z 2 )^ 



c 



v h,c(y ,Vi ,y 2 , x i x L,p) , (20) 

x (1 - e - iXLP+y z)(l - e - iXLP+{v ~-^ ] ) , (21) 
where the momentum fractions xl is defined as 

I 2 

°° L = 2p+q-z(l - z) ' (22) 



Note that the function Is t c{y~, Vi, 2/2", x, Xl, p) contains only phase factors. One can 
combine these phase factors with the matrix elements of the quark fields to define a 
special two-quark correlation function 



T, 



A(5,C) 



p + dy _ _ - 7 + _ - _ 7 + _ 
J =J 2n d Vi dy 2 (A\il> q (0)—ij} 9 (y )ip q (y 1 )^q{y 2 )\A) 



(23) 
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The contribution from quark-antiquark annihilation in the central-cut diagram in Fig. 4 
to the hadronic tensor can then be expressed as 




dz h 




C 2 F 2(l + z 2 ) A(5 , c) 
N c z(l-z) qq 



(x,x L ). 



(24) 



Contributions from all quark-quark (antiquark) double scattering processes can be cast 
in the above factorized form. 

The structure of the phase factors in I 5tC (y~, V\ , y^-, x, xl,P) is exactly the same as for 
gluon bremsstrahlung induced by quark-gluon scattering as studied in Ref. [18,19]. It 
resembles the cross section of dipole scattering and represents contributions from two 
different processes and their interferences. It contains essentially four terms, 



The first term corresponds to the so-called hard-soft processes where the gluon emission 
is induced by the hard scattering between the virtual photon 7* and the initial quark 
with momentum {x + xi)p. The quark then becomes on-shell before it annihilates with 
a soft antiquark from the nucleus that carries zero momentum and converts into a real 
gluon in the final state. The second term corresponds to a process in which the initial 
quark with momentum xp is on-shell after the first hard 7*-quark scattering. It then 
annihilates with another antiquark and produces two final gluons in the final state. In 
this process, the antiquark carries finite (hard) momentum XlP- Therefore one often refers 
to this process as double-hard scattering as compared to the first process in which the 
antiquark carries zero momentum. Set aside the change of flavors in the initial and final 
states, the double-hard scattering corresponds essentially to two-parton elastic scattering 
with finite momentum and energy transfer. This is in contrast to the hard-soft scattering 
which is essentially the final state radiation of the 7*-quark scattering and the total 
energy and momentum of the two final state gluons all come from the initial quark. The 
corresponding matrix elements of the two-quark correlation functions from these first two 
terms are called 'diagonal' elements. 

The third and fourth terms with negative signs in Is t c{y~ ,Vi ,2/2 ,x,Xl,p) are interfer- 
ences between hard-soft and double hard processes. The corresponding matrix elements 
are called 'off-diagonal'. The cancellation between the two diagonal and off-diagonal 
terms essentially gives rise to the destructive interference which is very similar to the 
Landau-Pomeranchuk-Migdal (LPM) interference in gluon bremsstrahlung induced by 
quark-gluon double scattering [18,19]. One can similarly define the formation time of the 




(25) 
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parton (quark or gluon) emission as 

T f = — —f ■ (26) 

In the limit of collinear emission (xl — > 0) or when the formation time of the parton 
emission, 77, is much larger than the nuclear size, the effective matrix element vanishes 
because 

h,c{v~,yx,V2,x,x L ,p)\ XL=Q ->b , (27) 



when the hard-soft and double hard processes have complete destructive interference. 

We should note that in the central-cut diagram of Fig. 4, the final state partons are two 
gluons. Therefore, in Eq. (20) the gluon fragmentation function in vacuum D g ^h{zh/ ' z) 
enters. If the other gluon (close to the 7*-quark interaction) fragments, the contribution 
to the semi-inclusive hadronic tensor is similar except that the corresponding effective 
"splitting function" should be replaced by 

z 1 — z) z 1 — z) 



As we will show in Appendix A-l, the two gluons in the quark-antiquark annihilation 
processes (central-cut diagrams) are symmetric when contributions from all possible an- 
nihilation processes and their interferences are summed. Therefore, one can simply mul- 
tiply the final results by a factor of 2 to take into account the hadronization of the second 
final-state gluon. 

In addition to the central-cut diagram, one should also take into account the asymmetrical- 
cut diagrams in Fig. 4, which represent interference between gluon emission from single 
and triple scattering. The hard partonic parts are mainly the same as for the central-cut 
diagram. The only differences are in the phase factors and the fragmentation functions 
since the fragmenting parton can be the final-state quark or gluon. These hard parts can 
be calculated following a similar procedure and one gets, 



a 2 x B fd£ 2 T rdz-. ,z h .C 2 F 2{l + z 2 ) 



-jjD , _ _ _ , a;x B [at? f dz z h 



Q 2 J l z T J z q y z' N c z(l-z) 

x Is,l(R) (V~ ' Vi ' V2 ■> x > x l , P) , (29) 
hAy~, VT, V2, x L ,p) = -e^^v- (1 _ e-W-ftr-tft-)) 

xe(yT-y;)9(y--yi), (30) 

hAV, VT, Th,x, x L ,p) = -e^+*^ + y- (1 - e'™^) 

x0(-y;)9(y;-y;). (31) 
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In the asymmetrical cut diagrams, the above contributions come from the fragmentation 
of the final-state quark. Therefore, quark fragmentation function D q ^h(zh/ z) enters this 
contribution. For gluon fragmentation into the observed hadron in this asymmetrical-cut 
diagrams, the contribution can be obtained by simply replacing the quark fragmentation 
function by the gluon fragmentation function D g ->h(zh/z) and replacing z by 1 — z. 
Summing the contributions from three different cut diagrams of Fig. 4, we can observe 
further examples of mixing (or conversion) of quark and gluon fragmentation functions. 
This medium-induced mixing was first observed by Wang and Guo [18] and is a unique 
feature of quark-quark (antiquark) double scattering among all multiple parton scattering 
processes. 

With the same procedure we can calculate contributions from all other cut diagrams 
of quark-quark (antiquark) double scattering at order G(a 2 s ), which are listed in Ap- 
pendix A-l. There are three types of processes: two annihilation processes, qq — > gg 
(central-cut diagrams in Figs. 5, 6, 7, 8 and 9), qq — »• q^ (central-cut diagram in Fig. 10) 
and quark-quark (antiquark) scattering, qqi(q~i) —> qqi(q~i) (central-cut diagram in Fig. 11). 
One also has to consider the interference of s and t-channel amplitude for annihilation 
into an identical quark pair, qq — > qq (central-cut diagrams in Figs. 12 and 13) and the 
interference between t and u channels of identical quark scattering qq — > qq (central-cut 
diagram in Fig. 14). 

Contributions from left and right-cut diagrams correspond to interference between the 
amplitude of gluon radiation from single 7*-quark scattering and triple quark scattering. 
The amplitudes of gluon radiation via triple quark scattering essentially come from ra- 
diative corrections to the left and right-cut diagrams of the lowest-order quark- antiquark 
annihilation in Fig. 3 (as shown in left and right-cut diagrams in Figs. 5, 6, 8, 9, 12, 
13, 15 and 16). Two other triple quark scatterings with gluon radiation, shown as the 
left and right-cut diagrams in Figs. 11 and 14, correspond to the case where one of the 
final state quarks, after quark-quark scattering, annihilates with another antiquark and 
converts into a final state gluon. 



4 Modified Fragmentation Functions 

In order to simplify the contributions from quark-quark (antiquark) scattering (annihi- 
lation), one can first organize the results of the hard parts in terms of contributions from 
central, left or right-cut diagrams, which are associated by contour integrals with specific 
products of ^-functions, 




Hg0(-y;)0(y- - „r) - #f %r - V^iV ~ Vi) 




(32) 
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These ^-functions provide a space-time ordering of the parton correlation and will restrict 
the integration range along the light-cone. For contributions from central, left and right- 
cut diagrams that have identical hard partonic parts, Hq^ c ' = = H^ c \ they will 
have a common combination of ^-functions that produces a path-ordered integral, 



y~ Vi 

dVi dyz=- I dy^dy^ e(-y^)6(y~ - ) - e(-y^)6(y^ - ) 



o o 



-6(y- - y^)9(y^ - y 2 ") (33) 



that is limited only by the spatial-spread y~ of the first parton along the light-cone 
coordinate. For a high-energy parton that carries momentum fraction xp + , y~ ~ l/xp + 
should be very small. Those contributions that are proportional to the above path-ordered 
integral are referred to as contact contributions (or contact interactions). 

Similarly, y{ — y^ is the spatial spread of the second parton and can only be limited 
by the spatial size of its host nucleon even for small value of momentum fraction. The 
spatial position of its host nucleon, y{~ + y^ , however, can be anywhere within the nu- 
cleus. Therefore, any contributions from double parton scattering that have unrestricted 
integration over y± and y^ should be proportional to the nuclear size of the target A 1 / 3 
and therefore are nuclear enhanced. In this paper, we will only keep the nuclear enhanced 
contributions and neglect the contact contributions. This will greatly simplify the final 
results for double parton scattering. 



4-1 qq — > g annihilation 



For the lowest order of quark- ant iquark annihilation in Eqs. (14)-(16), the hard parts 
from the three cut diagrams are almost the same except for the parton fragmentation 
functions. The central-cut diagram is proportional to the gluon fragmentation function 
while the left and right-cut diagrams are proportional to quark fragmentation functions. 
Rearranging the contributions from the three cut diagrams and neglecting the contact 
term that is proportional to the path-ordered integral as in Eq. (33), the total contribution 
can be written as 



E / dxT q f H \x,0)^2C F ^H^(x,p,q) 



dz h gJ N c ~ r Q* 

x [D g ^ h (z h ) - D q ^ h (z h )\ . (34) 

According to our definition in Eq. (13) of the twist-four correction to the quark fragmen- 
tation functions, the modification to the quark fragmentation function from the lowest 
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order quark-antiquark annihilation is then, 

^D^\z h ) = ^C F ^ 2 [D g _ h (z h ) - D^ h { Zh) ] T ^ ( ( ^ 0) • (35) 

Here the effective quark-antiquark correlation function T^ H ^ (x, 0) is defined as, 



T q ^ H \x,x L ) = J P^dy^dy^e^ + y--^ + ^^9(-y^9(y~ - yf) 

x (A\^ q {0)^ q {y-)^{yi)^Mv2)\A) , (36) 

with the antiquark carrying momentum fraction Xl. This two-parton correlation func- 
tion is always associated with double-hard rescattering processes. Similarly, we define 
three other quark-antiquark correlation matrix elements 

p + dy~ 

Vi 



x e{- y ;){A\$M\^y-)^M)\My;)\A) , (37) 
T q f~ L \x, x L ) = / t^ dyTdy -s*+^ + y-- i *^y--y^6{ y - - 

J Ztt 

x e(-yj)<A|^(o)^(y-)^«(yr)Y^*^2")l-A> , (38) 



T q f- R \x : x L ) = J P^-dy^dy^e^ + ^ + y—^6(y- 



Vx 



x 6(-y;)(A\^ q (0)^ q (y-)^M) , (39) 

that are associated with hard-soft rescattering and interference between double hard 
and hard-soft rescattering. In the first parton correlation (x, Xl), the antiquark q^ 

carries momentum fraction xl while the initial quark has the momentum fraction x. The 
two-parton correlation T^f\x, xl) corresponds to the case when the leading quark has 
x + Xl but the antiquark carries zero momentum. The two interference matrix elements 
are approximately the same for small value of xl and will be denoted as t£^(x, Xl)- 



4-2 qq —>■ qiqi annihilation 



Contributions from the next-to-leading order quark-antiquark annihilation or quark- 
quark (antiquark) scattering are more complicated since they involve many real and vir- 
tual corrections. The simplest real correction comes from qq — > q^ annihilation ^ q) 
[Fig. 10 and Eqs. (A-25) and (A-26)] which has only a central-cut diagram, 
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Cf o? s xb 
X E [D« 



dz 



z 2 + (1 - z?\ 



*h{z h /z) + D ?i _ h (z fc /z)] 



T nn (x ? x j 



(40) 



This kind of gg annihilation is truly a hard processes and thus requires the second an- 
tiquark to carry finite initial momentum fraction Xl- Furthermore, there are no other 
interfering processes. 



4-3 qqMi) -> qqMi) scattering 



Contributions from non-identical quark-quark scattering qq~i — > qqi (qi ^ q) are a little 
complicated because they involve all three cut diagrams (central, left and right) [Eqs. (A- 
28)-(A-32)]. One can factor out the ^-functions in the hard parts according to Eq. (32) 
and re-organize the phase factors in each cut diagram, 



j n c — e i ( x + x L)p + y~ m _ e -ix L p + V2 \n _ e -ix L p + {y~-yi)\ 

_ e i(x+x L )p+y- ^ _ e ~ix L p+y~ _ e ~ix L p+(y- -y~) _|_ e ~ix L p+ (y~ +y~ -y~ ) j . 
J n L — e i{x+x L )p+y~ M _ e -ix L p+(y--y~)^ 

_ e i(x+x L )p+y- ^ _ e -ix L p+y~ _ e ~ix L p+(y- -y~) _|_ e ~ix L p+y^ j . 

J u R = e i ( x + x L)P + V' (1 _ e ~iXLP + V2 \ 

_ e i(x+x L )p+y- ^ _ e -ix L p+y~ _ e ~ix L p+(y- -y~) _|_ e ~ix L p+ (y~ -yf )j 

such that the first three terms in each amplitude are the same. These three common phase 
factors will give rise to a contact contribution for all similar hard parts from the three 
cut diagrams, which we will neglect since they are not nuclear enhanced. The remaining 
part will have the following phase factors, 

7 n = e i( x + x L) P + y- [9(- y -)9{y- - y-) e - iXL P + ^ + y^-y^ 

-Ofa - yz)0(y- - yr)e- te " + *"" - - y-)e^^^r^\ . (42) 

Note that the phase factors of the last two terms in the above equation give identi- 
cal contributions to the matrix elements when intergated over and y^ as they dif- 
fer only by the substitution y^ <-> y\ — y~ ■ One therefore can combine them with 
9(— y 2)0(1/^ — yi)e~ lXLP+ ( y ~ :y i ) to form another contact contribution (path-ordered) 
which can be neglected. The final effective phase factor is then 

7 n = e ix P + y~ ~ ix LP + {y 2 -y 1)^1 _ e ix hP + V2^ _ ^43^ 
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Using the above effective phase factor, one can obtain the effective modification to the 
quark fragmentation function due to quark-antiquark scattering, qqi — > qq i: 



AD 



q—*h 



Cf o? s x b f d£j< r dz 

^■qi J Z 



n c g 2 

+ Dg^ h (z h /Z] 



E 

i 

1 + (1-^) 2 



D q ^ h (z h /z) 



1 + z* 



T, 



{l-zf 
f x,x L ) 



<!>! 



+ [D^ h (z h /z)) - D g ^ h {z h /z)\ 



2 rpA(HS), 



l + (l- Z ) 2 T^>(x,X L ) 



n c g 2 

+ Dg % ^ h {z h /z) 



Cf o? s xb f di\ r dz ^ j 

a J 

l + (l-z) 2 



D q ^ h (z h /z) 



1 + z 2 



(l-,) 2 

)f X,X L ) 



qq 



+ [D rH ^ h (z h /z) - D g ^ h (z h /z)} 



f q A (x) 

2r r A{SI), 



where three types of two-parton correlations are defined: 



(44) 



Tqqi \ x i x L)—T qq ^ \x,Xl) T qq -.\x,XL) 

= f P+dy dy^dy^e ixp+v ~- iXLP+ ^" y ^\l - e iXLP+y *) 
J 2ir 

x(il|^(0)^(y-)^(yr)^«(y 2 -)|i4>e(-ife)fl(y- - yf) 



Tq q \ (x, Xl) — T qq . (x,Xl) T qq . (x,Xl) 



p + dy~ 
2tt 



-dy 1 dy 2 e 



i(x+x L )p+y 



(1-e 



■ix L p^y 2 



x<A|^(o)^(r)^*(yr)^*(i/2")l> 4 > e (-%) e (i/" -vt) > 



r A(HS) 

l qqi ' 



x,x L ) = T qq {HI \x, x L ) + T qq f n (x, x L ) 



Msi) i 



p + dy~ 
2vr 



■dy^dy^e i(x+XL)p+y (1 - e~ iXLP+y ^ ) 



x (i _ e ^LV+{y--yn)e(-y-)6( : y- - yf) 



x(A|^(o)l-^(y-)4(2/r)Y^^)l^- 



,7 



(45) 



(46) 



(47) 



One can similarly obtain the modification of quark fragmentation from non-identical 
quark-quark scattering by replacing q\ — > qi in Eq. (44), 
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Cf a'jxB f d&r t dz 

w 



n c q 2 

+ D qi ^ h (z h /z) 



E 

<n+<i 
l + (l-z) 2 ' 



D q ^ h (z h /z) 



1 + z 2 

(T^zf 



+ [Dq.^ h (z h /z) - D g ^ h (z h /z)} 



l + (l-zfT q A q f\ Xl x L ) 
z 2 fMx) 



(4* 



The two-quark correlations, T^ HI \x, x£) and T^ SI \x, x£) can be obtained from T qq . (x,xl) 
and Tggf^ (x,x L ), respectively, by making the replacements ifj qi (y 2 ) -> $ qi (y 2 ) and $ q . {y x ) -> 
^(yi) in Eqs. (45) and (46), 



T^ HI) (x,x L ) = J ^-^-dy x dy 2 e ixp+y - iXLP+ ^-^\l - e iXLP+y *) 

x(A|^(0)^ fl (y-)^(yJ)^ tt (yr)|A)fl(-i/ a -)e(i/- - yf) , (49) 
T# 5i WiO = / ^^dy^dy^+^y- (1 - e - ix ^) 

x^l^^^^y-^fy^Y^^I^^-^^-yr). (so) 

and T# flS >(z, x L ) = T^(x, z £ ) + T^)(x, x £ ). 

Note that the contribution from fragmentation of quark q^ or antiquark q~i only comes from 
the central-cut diagram. This contribution is positive and is proportional to T^) (x, xl) + 
T^ SI \x,Xl), containing all four terms: hard-soft, double-hard and both interference 
terms . The gluon fragmentation comes only from the single-triple interferences (left and 
right-cut diagrams). Its contribution is therefore negative and partially cancels the pro- 
duction of qi(<ji) from the hard-soft rescattering. The cancellation is not complete since 
the gluon and quark fragmentation functions are different. The structure of this hard-soft 
rescattering (quark plus gluon) is very similar to the lowest order result of qq — »• g in 
Eq. (35). It contributes to the modification of the effective fragmentation function but 
does not contribute to the energy loss. The energy loss of the leading quark comes only 
from double-hard rescattering, since the leading quark fragmentation comes both from 
the central-cut and single-triple interferences, and the single-triple interference terms can- 
cel the effect of hard-soft scattering for the leading fragmentation. Its net contribution is 
therefore proportional to T^S^JJ . Since the double-hard rescattering amounts to elastic 
qqi(q~i) scattering, the effective energy loss is essentially elastic energy loss as shown in 
Ref. [34]. There is, however, LPM suppression due to partial cancellation by single-triple 
interference contributions. For long formation time, l/xip + ^> Ra, the cancellation is 
complete. Therefore, LPM interference effectively imposes the lower limit xl > 1/p + Ra 
on the fractional momentum carried by the second quark (antiquark). 
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4-4 QQ ~~ *■ QQ scattering 



For identical quark-quark scattering, qq — > qq, one has to include both t and ^-channels, 
their interferences, and the related single-triple interference contributions. Using the same 
technique to identify and neglect the contact contributions, one can find the correspond- 
ing modification to the quark fragmentation function from Eqs. (48) and (A-45)-(A-49), 



^™\z h 



C F a 2 s x B r df T } dz / T^ HS \x, x L 



rp J Z 



x [Dg^ h (z h /z)) - D g ^ h {z h /z)} 

+ 



:i-z) 



n ( I \ ( 1 + z * 1 1 



+ Dg^ h (z h /Z) 



:i-z) 



Nrz(l-z) 



N r z(l - z) 



T^ HI \x,x L ) 



C F o? s x B rdf T }dz\ T^{x,x L ) 

N c Q 2 J t\ J z \ fMx) ^-m^h^h 

z h 



' z h/z) 



-D g _+ h (z h / z)\ +D q ^ h (z h / z)P qq ^ qq (z) 



T q A q ( HI \x,x L ] 



(51) 



where the effective splitting functions are defined as 



P {s) (z) 

qq^qqx J 



1 + (1 - zf 



P 



qq-^qq 



>) 



:i-z) 



+ 



i i 

l + z 2 



[l-zf N c z(l-z) 



(52) 
(53) 



4-5 qq — > qq, gg annihilation 



The most complicated twist-four processes involving four quark field operators are quark- 
antiquark annihilation into two gluons or an identical quark-antiquark pair. We have to 
consider them together since they have similar single-triple interference processes and 
they involve the same kind of quark-antiquark correlation matrix elements, T qq ^(x, Xl), 
(i = HI, SI,HS). 

For notation purpose, we first factor out the common factor (C F / N c )a 2 x b / Q 2 / f q (x) and 
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the integration over £ F and z and define 

C F o?xb 



AD^^\z h 



N c Q 2 L A (x) 



fj^AD^ 99 ^(z h ,z,x,x L ). 

T Zh 



(54) 



After rearranging the phase factors and identifying (by combining central, left and right 
cut diagrams) and neglecting contact contributions we can list in the following the twist- 
four corrections to the quark fragmentation from the hard partonic parts of each cut 
diagram (see Appendix A): 

Fig. 5 (t-channel qq — > gg), 



l + (l-z) 



1 + z 2 



z{l — z) z(l — z 
1 + z 2 



A{HI) 



qq 



D g ^ h (z h /z)2C F 
+ [D g ^ h (z h /z) - D q ^ h (z h /z)} 2C F T^ {SI) (x, x L 
Fig. 6 (interference between u and t-channel of qq — > gg) 



x,x L ) 



(55) 



AD 



- qq 



(qq^ 9 g,qq)_ n , / x ~4 (C F ~ C A /2) A(HI) 

q^h(6) -Ug^ h {Z h /Z) — - - ^ 

+ [D g ^ h (z h /z) - D q ^ h (z h /z)} 



-2{C f — C a /2) a{sj) ( n. 

1 qq \X, X F ) , 



zl-z) 



Fig. 7 (s-channel of qq — > gg), 



AD 



= D g ^ h {z h /z)AC A ^ ~ ' T q T'(x,x L ) ; 



1 - z + z 2 ) 2 ^ A[H) 



Figs. 8 and 9 (interference of s and t-channel qq — > gg), 



(56) 



(57) 



= D^ h (z h /z)(-2C A ) 



xT q f HI \x,x L ) + C A 



1 + z 3 l + (l-z) 



-D g ^ h (z h /z 



1 + (1 ~ *) 
z(l-z) 



z(l-z) 
D q ^ h (z h /z) 

,A(I2) 



z{l-z) 

1 + Z 3 
2(1 - *) 



x [T^\x,x L )-T^>(x,x L )] 



(58) 



Fig. 10 (s-channel of qq — > qq), 
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= [D q ^ h (z h /z) + D^ h (z h /z)] [z 2 + (1 - z) 2 } 



q^/i(10) 



xT$ H \x,x L ) 



Fig. 11 (i-channel of qq 



similar to Eq. (44), 



(59) 



AD 



(qq^gg,qq) _ n ( i s - . ~ A yni } , , 
q^h(\\) -Uq^hKZh/Z) , 1 _ z j 2 *qq \X,Xl) 



1 + Z 2 A(HI) 



+Dg^ h (z h /z) 

-Dg^ h (Z h /z) 



l + (l-z) 2 A{HS) 



qq 



X,X L , 



:i-z) 



-,A(SI). 



- qq 



x,x L ) 



(60) 



Figs. 12 and 13 (interference between s and t-channel qq — > qq), 



a n{qq-^99,qq) 

^-^9^/1(12+13) 



-4(C F - CU/2) 



D q ^ h (z h /z) — 



+ Dg_ h (Z h /z] 



(1-zf 



Tqq {x, Xl) 



-2(C F - C A /2) 



D q ^ h (z h /z)- h D g ^ h (z h /z) 



(1-zf 



1 - z 



X \Tqq {I2 \x,X L ) ~ T qq {I \x,X L )\ ] 

Figs. 15 and 16 (two additional single-triple interference diagrams), 



(61) 



^%^h(15+16) 



-2C f 



D q ^ h (z h /z) 



1 + z 2 



1-2 



+ Dg^ h (Z h /z) 



l + {l-zf 



,A(/2) 
qq \ 



(62) 



Most processes involve both T a ^ hi \x,xl) for double-hard rescattering with interference 
and T A ( SI \x,x L ) for hard-soft rescattering with interference. All the s-channel (Figs. 7 
and 10) processes involve double-hard scattering only. Therefore, they depend only on the 



T, 



<i<i 



X,X L) 



nA(HI) , 

■qq 



x,xl) -\-T q q{x,xi). For interference between single and triple 
scattering (left and right-cut diagrams in Figs. 8, 9, 12 13, 15 and 16), where a hard 
rescattering with the second quark (antiquark) follows a soft rescattering with the third 
antiquark (quark), only interference matrix elements, T qq (x, Xl) and T qq \x, Xl), are 
involved. Here, 



Tqq (%, % l) 



p + dy~ 
2vr 



•dy 1 dy 2 e 
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p + dy~ 
2vr 



■dy{dy^e ixp+y +ix ^ + ^ -Vi) 



x(A|^(o)^(y-)^(yr)^(y 2 -)i^)^(-^)^(r-yr 



(63) 



is a new type of interference matrix elements that is only associated with this type of 
single-triple interference processes. One can categorize the above contributions according 
to the associated two-quark correlation matrix elements and rewrite the above contribu- 
tions as, 



= T ^ HI) (*» x ^ [D g ^ h {z h /z)P q ^ gg {z) + D q ^ h (z h /z)P qq -^ qq ( 



+D q ^ h (z h / z)P qq ^ qq (l ~ z)} 



^^q^hisi) — 1 qq \ X i X L) 



a 



z(l - z) 
xD g ^ h (z h /z) - D q ^ h (z h /z 



2C r 



i + (i- z y 



1 - z 

c A 



+ D q ^ h (z h /z) 



z{l-z) 
l + (l-z) 2 



+ 2C t 



1 - z 



2C f 



r A(SI) 
L qq 1 

1 + z 2 



z(l-z) 

- T q q (x, Xl) 



x,x L ) {[D q ^ h (z h /z) - D g ^ h (z h /z)} [P£l qq (z) 

+ [D q ^ h (z h /z) - D q ^ h (z h /z)} 1 + (1 ~ Z)2 

z 



^jjqq^qq,99 _ rpA(i) 



r 4(l-z + z 2 ) 2 -l (1-z) 2 

^ A JZ \ ZLyp 



z(l — z) 

x Dg-+ h (z h /z) + [z 2 + (1 - zf]D q ^ h {z h /z)\ 



+ D q ^ h (z h /z) 



z 2 + (1 - z) 2 - 



C 



2C f 



+ T qq {I2 \x 1 x L ) i.D q ^ h (z h /z) 
+ D g ^ h (z h /z) 



z{l-z) l \-z 
C a 2C p 



+ 



z{l-z) 1-z 
C a 2C p 
z(l — z) z 



where P qq l^ qq (z) is given in Eq. (52) and the effective splitting functions for qq 
qq — > qq are defined as 



(64) 



(65) 



(66) 
gg and 



Pqq-^gg( z ) 



Pqq^qq(z) 



z r _ 2Ca[z2 + (1 _ z)2] . 



z(l-z) 
(l-^) 2 4 



1 + z 2 



N c l 



(67) 
(68) 
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which come from the complete matrix elements of qq — > gg and qq — > qq scattering (see 
Appendix A-3). Again, double-hard rescattering corresponds to the elastic scattering of 
the leading quark with another antiquark in the medium and the interference contribu- 
tions. The structure of the hard-soft rescattering contribution we identify above shows 
the same kind of gluon-quark (or quark-antiquark) mixing in the fragmentation functions 
and does not contribute to the energy loss of the leading quark. The unique contributions 
in the qq — > qq, gg processes are the interference-only contributions. They mainly come 
from single-triple interference processes in the multiple parton scattering. 



5 Modification due to quark-gluon mixing 



We have so far cast the modification of the quark fragmentation function due to quark- 
quark (antiquark) scattering (or annihilation) in a form similar to the DGLAP evolution 
equation in vacuum. In fact, one can also view the evolution of fragmentation functions in 
vacuum as modification due to final-state gluon radiation. In both cases, the modification 
at large is mainly determined by the singular behavior of the splitting functions for 
z — > 1, whereas the modifications at mall Zh is dominated by the singular behavior of the 
splitting function for z — > 0. 

Let us first focus on the modification at large Zh- A careful examination of the contribu- 
tions from all possible processes shows that the dominant modification to the effective 
quark fragmentation function comes from the t-channel of double hard quark-quark scat- 
tering processes, 



Cp a 2 x B „ 



N c Q 2 



di\ j dz^^ ,Zh. 

£rri J Z Z 



T^ HI \x,x L ] 



<^L+6(l-z)A gi (f T ) 



z)t 



df 



T " 



dz 



C F 2 ^ r ar T r az z h 
Jf a s L / — / — D ^h{— 



z 



+ S(l - z)A qi (f T ) 



(69) 



where the summation is over all possible quark and antiquark flavors including q% = q,q 
and A qi [ij) represents the contribution from virtual corrections. We have expressed the 
modification in a form that it is proportional to the matrix elements xlT^ hi \x, xl) / f t 



A 1 / 3 xl/ ( ^(xl) as compared to the modification from quark-gluon scattering where the 
corresponding matrix element [Eq. (9)] is T^ HI \x, xl) / f^(x] 



) ~ A l ^x L G N {x L ). Here, 
fq (x) and G N (x) are quark and gluon distributions, respectively, in a nucleon. This 
leading contribution to the modification from quark-quark scattering is very similar in 
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form to that from quark-gluon scattering [see Eq. (7)]. However, it is smaller due to the 
different color factors Cf/Ca = 4/9 and the different quark and gluon distributions, 
fq i {xi) and G n (xl) in a nucleon. Because of LPM intereference, small angle scattering 
with long formation time ry = 1/xlp + is suppressed, leading to a minimum value of 
xl > xa = ^/tunRa = 0.043 for a Kr target. For this value of xl, the ratio 

V f N (xr O 2 ) 

%%^P> 1-40/1.85 - 0.75, (70) 



at Q 2 = 2 GeV 2 according the CTEG4HJ parameterization [35]. Therefore, one has to 
include the effect of quark-quark scattering for a complete calculation of the total quark 
energy loss and medium modification of quark fragmentation functions. 

In a weakly coupled and fully equilibrated quark-gluon plasma, quark to gluon number 
density ratio is p q j p g = rif(3/2)N c / (N 2 — 1) = 9n//16. An asymptotically energetic 
jet in an infinitely large medium actually probes the small x = (q^)/2ET regime, where 
quark- ant iquark pairs and gluons are predominantly generated by thermal gluons through 
pQCD evolution. In this ideal scenario one expects N q /N g ~ 1/4Ca = 1/12 and therefore 
can neglect quark-quark scattering. The modification of quark fragmentation function will 
be dominated by quark-gluon rescattering. However, for moderate jet energy E « 20 GeV 
and a finite medium L ~ 5 fm, parton distributions in a quark-gluon plamsa are close to 
the thermal distribution. In particular, if quark and gluon production is dominated by 
non-perturbative pair production from strong color fields in the initial stage of heavy-ion 
collisions [36], the quark to gluon ratio is comparable to the equilibrium value. In this 
case, we should take into account the medium modification of the quark fragmentation 
functions by quark-quark scattering. 

An important double hard process in quark-quark (antiquark) scattering is qq — > gg 
[Eq. (64)]. In this process, the annihilation converts the initial quark into two final gluons 
that subsequently fragment into hadrons. This will lead to suppression of the leading 
hadrons not only because of energy loss (energy carried away by the other gluon) but 
also due to the softer behavior of gluon fragmentation functions at large z^. Even though 
the leading behavior of the effective splitting function [Eq. (67)] 



z 2 



P q ^ gg {z) « 2C F ±_ ' (71) 



is not as dominating as that of t-channel quark-quark scattering, it is enhanced by a color 
factor 2Cf = 8/3. One expects this to make a significant contribution to the medium 
modification at intermediate Zh- 

In high-energy heavy-ion collisions, the ratios of initial production rates for valence 
quarks, gluons and antiquarks vary with the transverse momentum p?- Gluon production 
rate dominates at low p T while the fraction of valence quark jets increases at large px- 
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Quarks are more likely to fragment into protons than antiprotons, while gluons fragment 
into protons and antiprotons with equal probabilities. Therefore, the ratio of large pt 
antiproton and proton yields in p + p collisions is smaller than 1 and decreases with px 
as the fraction of valence quark jets increases. Since gluons are expected to lose more 
energy than quark jets, one would naively expect to see the antiproton to proton ratio 
p/p becomes smaller due to jet quenching. However, if the quark-gluon conversion due 
to qq — > gg becomes important, one would expect that the fractions of quark and gluon 
jets are modified toward their equilibrium values. The final p/p ratio could be larger than 
or comparable to that in p + p collisions. Such a scenario of quark-gluon conversion was 
recently considered in Ref. [37] via a master rate equation. 

The mixing between quark and gluon jets also happens at the lowest order of quark- 
antiquark annihilation as shown in Fig. 3. At NLO, all hard-soft quark-quark (antiquark) 
scattering processes have this kind of mixing between quark and gluon fragmentation 
functions. Their contributions generally have the form, 



C 



F OiZXB 



E 



di\ 



dz 

z 



~h 



xP, 



[) fftx) 



(72) 



where again the summation over the quark flavor includes = q, q. This mixing does 
not occur on the probability but rather on on the amplitude level since it involves in- 
terferences between single and triple scattering. Therefore, this contribution depends on 
the difference between gluon and quark fragmentation functions [Eq. (35)] and can be 
positive or negative in different region of Zh- Nevertheless, they contribute to the mod- 
ification of the effective quark fragmentation function and the flavor dependence of the 
final hadron spectra. 



6 Flavor dependence of the medium modified fragmentation 



Summing all contributions to quark-quark (antiquark) double scattering as listed in Sec- 
tion 4, we can express the total twist-four correction up to 0{a 2 s ) to the quark fragmen- 
tation function as 



AD q ^ h (z h ) 



N r . 



Q 2 
i 



2 [Dg^ h (zh) - Dq^ h (Zh)] 



x,0) 



fa 



^/f 1 j d iY,n^ h {z h /z)p«U^ qa 



T g A ^(x,x L ) 



a,b,i 



(73) 
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where the summation is over all possible q + a — > b + X processes and all different matrix 
elements T^ l '(x, xl) (i = HI, SI, 1, 12), which will be four basic matrix elements we will 

use. The effective splitting functions Pq}^ b {z) are listed in Appendix A-2. One should also 
include virtual corrections which can be constructed from the real corrections through 
unitarity constraints [18]. 

Similarly, we can also write down the twist-four corrections to antiquark fragmentation 
in a nuclear medium, 



An ( \ Cf o asXB lorn ( \ n t M T ti °) 
AD^ h (z h ) = j^^-qT { 2 [Dg^ h (z h ) - D^ h (z h )\ fA ^ 



where the matrix elements T^ 1 \x,Xl) and the effective splitting functions P^_^ b (z) can 
be obtained from the corresponding ones for quarks. Given a model for the two-quark 
correlation functions, one will be able to use the above expressions to numerically evaluate 
twist-four corrections to the quark (antiquark) fragmentation functions. In this paper, 
we will instead give a qualitative estimate of the flavor dependence of the correction in 
DIS off a large nucleus. 

For the purpose of a qualitative estimate, one can assume that all the twist-four two-quark 
correlation functions can be factorized, as has been done in Refs. [18,19,23,33], 



' dy7dy^e^ + y~ + ^ p+ ^-^6(-y^6(y- - y7 



J 2tt 

x<i4|^(o)^ 9 (y-)^(yr)Y^«(yi")l-A> 

« — f£M , (75) 

Xa 

J 2tt 

K—tftofSMe-*''*, (76) 
xa 

where xa = ^/m N RA, m N is the nucleon mass, Ra the nucleus size, f^(x 2 ) is the 
antiquark distribution in a nucleon and C is assumed to be a constant, parameterizing 
the strength of two-parton correlations inside a nucleus. The integration over the position 
of the antiquark (y7 + y~)/2 in the twist-four two-quark correlation matrix elements gives 
rise to the nuclear enhancement factor 1/xa = ttimRa = 0.21A 1 / 3 . 
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We should note that we set kx = for the collinear expansion. As a consequence, the 
secondary quark field in the twist-four parton matrix elements will carry zero momentum 
in the soft-hard process. Finite intrinsic transverse momentum leads to higher-twist cor- 
rections. If a subset of the higher-twist terms in the collinear expansion can be resummed 
to restore the phase factors such as e tXTP+y , where xt = (k^) /2p + q~ z(l — z), the soft 
quark fields in the parton matrix elements will carry a finite fractional momentum xt- 

Under such an assumption of factorization, one can obtain all the two-quark correlation 
matrix elements: 



T q ^ HI \x, x L ) « — tf(x) f»{x L + x T )[l - e-*ll*\ (77) 

Xa 

T^ SI \x,x L )*^tf(x + x L )f»(x T )[l-e-*l/*A), (78) 
xa 

T q f\x,x L )*T q f 2 \x,x L ) » —ft{x + x L )f»{xr)e-*.l** 

xa 

^—ffWf^XL + x^e-^. (79) 
xa 

In the last approximation, we have assumed x^ ~ xt <C x. Similarly, one can obtain 
T^\x,xl), Tq q f' \x,xl) and T^'{x,Xi). With these forms of two-quark correlation 
matrix elements, we can estimate the flavor dependence of the nuclear modification to 
the quark (antiquark) fragmentation functions. 

The lowest order corrections [C(a s )] are very simple 



AD^(z h ) ex CA^[D g ^ h {z h ) - D q ^ h (z h )}f q N (x T ) , (80) 



AD^(z h ) cx CA^[D g ^- h (z h ) - D^ h (z h )]f Q N (x T ) . (81) 



We consider the dominant contribution from the fragmentation of a quark (antiquark) 
which is one of the valence quarks (antiquarks) of the final particle h (antiparticle h). 
The gluon fragmentation functions into h and h are the same. For large Zh, the gluon 
fragmentation function is always softer than the valence quark (antiquark) fragmenta- 
tion [38]. Therefore, the lowest order twist-four corrections are always negative for large 
Zh, leading to a suppression of the valence quark (antiquark) fragmentation function, 
D qv ^h{ z h) [Dq v ->h{zh)]- Consider those quarks that are also valence quarks of a nucleon: 



n = udd 

p = uud ,p = uud , (82) 

K + = us ,K~ = us . (83) 

7r + ,7r°,7r~ = ud ,(uu —dd)/V2 ,du . (84) 

One can find the following flavor dependence of the lowest order twist-four corrections 
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to the quark (antiquark) fragmentation functions 



^Df v °} h {z h ) -\AD™(z h )\ 



> 1 , (85) 



or 



where R q L v °}h ls the corresponding leading order suppression of the fragmentation function 
at large Zh for proton (anti-proton) and K + (K~). Since pions contain both valence 
quark and antiquark, the suppression factors should be similar for all pions. For iy > 
0.043, u(x)/u(x) > 3 and d(x)/d(x) > 2 [35]. Therefore, the modification of antiquark 
fragmentation functions due to quark- antiquark annihilation is significantly larger than 
that of a quark. 



The flavor dependence of the NLO results are more complicated since they involve scatter- 
ing with both quarks and antiquarks in the medium. One can observe first that effective 
splitting functions (or quark-quark scattering cross section) are the same for the t-channel 
qq' — > qq' and qq' —> qq' (q' ^ q) scatterings, 



For identical quark-quark scattering or quark-antiquark annihilation, one can separate 
the qq annihilation splitting functions (or cross sections) into singlet and non-singlet 
contributions by singling out the t-channel contributions, 



P$U*) = P«U Z ) + Ai^(z), (88) 
P^ b {z) = P^ b {z) + AI^! b (z). (89) 

These singlet contributions to the modified fragmentation functions are, 



T b,q',i 

xV?(*t) + $(xt)]C® , (90) 

T b,q',i 

*\f?M + f?M]C®, (91) 
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where the summation over q' now includes q'—q and C^'(xl) are flavor-independent 
functions determined from Eqs. (77)-(79), 



c (Hi) =c (si) = c(x L )(l - e- x i /x *), 
C (i) = C (m = c{x L )e- x V x A , (92) 

and C{xl) is a common coefficient that is a function of x F . Using P^)_^{z) = PgJ^ b {z) 
, one can conclude that the singlet contributions to the modified quark and antiquark 
fragmentation functions are the same, ADg^ h LO \zh) = AD S -^^°\zh)- 

The non-singlet contributions, mainly from s-channel and s-t interferences, are, 



A<if °\z h ) oc |MV3 Db ^ h ® AP»%(z h )f»(x T )CU , (93) 

AD^f°\z h ) oc /< E ® AP»%(z h )f?(x T )CU , (94) 



" T b,i 



where again AP ^ b (z) = AP. _k{z) due to crossing symmetry. We have listed all non- 



vanishing nonsinglet splitting functions AP^^ b (z) in Appendix A-2. 

We again consider the limit z^ — > 1. In this region the convolution in the modified 
fragmentation function is dominated by the large z — > 1 behavior of the effective split- 
ting functions. From the listed AP^^ b (z) in Appendix A-2, we can obtain the leading 
contributions, 

i 1 Z 

C^AP^z) « 2 [2C F + C F {1 - e~^A) + C A e~ x2 ^] (95) 



where we have also neglected terms proportional to 1/N C . All AP^^q(z) are non-leading 
in the limit z — > 1 and therefore can be neglected. With these leading contributions, 
the non-singlet modification to the quark and antiquark fragmentation functions can be 
estimated as 



AD N q i N h LO \z h )^A^ J 



dz 



D 



g^h 



x (CV(1 - e-*M) + C A e< /a ^) + 6(1 - z)A 1 (£ T ) 



C{x L ) 



D 



g->/i 



D 
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X 



? 3 XL } +S(1-z)A 2 (£ t ) 



l-z 



+ 




(96) 



7 -7i 

x (C F {1 - e-^A) + C A e- x i /x A) + 5(1 - z)A 1 {£ T ) 
+ 

C{x L ) 



M? "Ml 



X 



2C f 



+ 5(l-z)A 2 (£ T ) 



where Ai(#r) and ^{It) are from virtual corrections, 



f q N M , 



(97) 



Al ^) = / iZT z { C f C ^l)\ z =i - [C F {1 - e-*V*A) 
o 

+ C A e-*V*A]C{x L )} 
A 2 (£ r ) = J ^fC F [C{x L )\ z=l - C{x L )\ . 



(98) 
(99) 



Because of momentum conservation, C(xl) = when xl oo for z = 1. Therefore, 
the above virtual corrections are always negative. At large Zh, these virtual corrections 
dominate over the real ones. 

There are two kinds of non-singlet contributions in the expressions given above. One that 
is proportional to gluon fragmentation functions is due to quark-antiquark annihilation 
into gluons which then fragment. The fragmenting gluon not only carries less energy than 
the initial quark but also has a softer fragmentation function, leading to suppression of 
the final leading hadrons. The second type of contributions is proportional to D g ^h(zh) — 
D q ^h{zh) and therefore mixes quark and gluon fragmentation functions, similarly as the 
lowest order quark-antiquark annihilation processes [see Eqs. (80) and (81)]. Since a gluon 
fragmentation function is softer than a quark one, the real corrections from this type of 
processes are positive for small Zh and negative for large Zh- The virtual corrections 
have just the opposite behavior. Therefore, the second type of contributions will reduce 
the total net modification. For intermediate values of Zh where 2D g ^ h (z h ) > D q ^ h (zh), 
the net effect is still the suppression of the effective fragmentation functions for leading 
hadrons. 



Since f^(xx) > /^(^t), we can conclude that the LO and NLO combined non-singlet 
suppression for antiquark fragmentation into valence hadrons is larger than that for quark 
fragmentation into valence hadrons. This qualitatively explains the flavor dependence of 
nuclear suppression of leading hadrons in DIS off heavy nuclear targets as measured by 
the HERMES experiment [25,26]. The ratio of differential semi-inclusive cross sections 
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for nucleus and deuteron targets were used to study the nuclear suppression of the frag- 
mentation functions. It was observed that suppression of leading anti-proton is stronger 
than for leading proton and K~ suppression is stronger than K + . In the valence quark 
fragmentation picture, the leading proton (K + ) is produced mainly from u, d (u) quark 
fragmentation while anti-protons come primarily from u, d (u) fragmentation. Therefore, 
HERMES data are consistent with stronger suppression of antiquark fragmentation. 

Since gluon bremsstrahlung and the singlet qqi(q~i) scattering also suppress quark and 
antiquark fragmentation, but independently of quark flavor, one has to include all the 
processes in order to have a complete and quantitative numerical evaluation of the flavor 
dependence of the nuclear modification of the quark fragmentation functions. Further- 
more, the NLO contributions are proportional to a s ln(Q 2 )/27r. They are as important 
as the lowest order correction for large values of Q 2 . In principle, one should resum these 
higher order corrections via solving a set of coupled DGLAP evolution equations, in- 
cluding medium modification for gluon fragmentation functions. The contributions from 
quark-quark (antiquark) scattering derived in this paper will be an important part of the 
complete dscription. Detailed numerical study of the effect of quark-quark (antiquark) 
scattering will be possible only after the completion of this complete description in the 
future. 



7 Summary 

Utilizing the generalized factorization framework for twist-four processes we have stud- 
ied the nuclear modification of quark and antiquark fragmentation functions (FF) due 
to quark-quark (antiquark) double scattering in dense nuclear matter up to order 0(a 2 ). 
We calculated and analyzed the complete set of all possible cut diagrams. The results 
can be categorized into contributions from double-hard, hard-soft processes and their 
interferences. The double-hard rescatterings correspond to elastic scattering of the lead- 
ing quark with another medium quark. It requires the second quark to carry a finite 
fractional momentum xjj. Therefore, the energy loss of the leading quark through such 
processes can be identified as elastic energy loss at order 0(a 2 ). The quark energy loss 
and modification of quark fragmentation functions are dominated by the t-channel of 
quark-quark (antiquark) scattering and are shown to be similar to that caused by quark- 
gluon scattering. The contribution from quark-quark scattering is smaller than that from 
quark-gluon scattering by a factor of Cf/Ca times the ratio of quark and gluon distribu- 
tion functions in the medium. We have shown that such contributions are not negligible 
for realistic kinematics and finite medium size. The soft-hard rescatterings mix gluon and 
quark scattering, in the same way as the lowest order qq — > g processes. Such processes 
modifies the final hadron spectra or effective fragmentation functions but do not con- 
tribute to energy loss of the leading quark. For qq — > qq, gg processes, there also exist 
pure interference contributions mainly coming from single-triple-scattering interference. 

With a simple model of a factorized two-quark correlation functions, we further investi- 
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gated the flavor dependence of the medium modified quark fragmentation functions in a 
large nucleus. We identified the flavor dependent part of the modification and find that 
the nuclear modification for an antiquark fragmentation into a valence hadron is larger 
than that of a quark. This offers an qualitative explanation for the flavor dependence of 
the leading hadron suppression in semi-inclusive DIS off nuclear targets as observed by 
the HERMES experiment [25,26]. 
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A-l Hard partonic parts for quark-quark double scattering 

In Section 3 we have discussed the calculation of the hard part of one example cut-diagram 
(Fig. 5) in detail. In this appendix we list the results for all possible real corrections to 
quark-quark (antiquark) double scattering in the next-to-leading order 0(a^). There are 
a total of 12 diagrams as illustrated in Figs. 5-16. For the purpose of abbreviation, we 
will suppress the variables in the notations of partonic hard parts 



We first consider all qq — > gg annihilation diagrams with different possible cuts. The 
contributions of Fig. 5 are: 



H =H (y ,y 1 ,y 2 ,x,p,q,z h ) , 



(A-l) 



and phase factor functions 



i = Hv ,yi ,V2,x,,xl,p) ■ 



(A-2) 




1 + (1-*) 2 1 ci 
z(l -z) \ N c ' 



(A-3) 
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Fig. 5. The t-channel of qq — ► gg annihilation diagram with three possible cuts, central(C), 
left(L) and right(R). 




i i i 

Fig. 6. The interference between t and ^-channel of qq — > gg annihilation. 

h,c =o(-y;)0(y--yi)e i{x+XL)p+y - 

x(l - e - iXLP+y z)(l - e - iXLP+{v '~^ ] ) , (A-4) 



n 5,L(i?) 



h,L ■ 



o> 2 s xb f f dz— 



5,L(R) 



D g ^ h (z h /z)2 



1 + z 2 



-D g ^ h (z h /z)2 



z(l-z) 
l + (l-z) 2 



-8(-y 2 )%2 - Vi ) e 



g-+ nK ~ nj ~,- z[l _ z) 

i(x+x L )p+y~ (-\ _ e -ix L p+(y 

i(x+XL)p + y 



Nr. 



(1 — e - tx LP y 2 



(A-5) 

(A-6) 
(A-7) 



Here we have included the fragmentation of both final-state partons. 



The contributions from Fig. 6 are: 
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Fig. 7. The s-channel of qq — > gg annihilation diagram with only a central-cut. 
i 



■rrD _atx B rdt% fdz 7 

-2 C F (C F -C A /2) 



'T 



z)z N c 



i 



— d ol\xb f d£^ f dz- 
H, 



Qi J l\ 



/dz— 
— h,L(R) [ D g^h(z h /z) + Dq-> h (z h /z)] 



(A-8) 



-2 C F (C F -C A /2) 

(l-z)z N c ' lAyj 
h,c = 0(-y^)e(y- - y-) e ^+^ + y- 

x (1 - e ~ lXLP+y ^){\ - e - iXLP+ ( y ~-^) , (A-10) 

h, L = -%r - y^iv - yT)e i{x+XL)p+y ~ (i - e-^y-y^) . (a-ii) 

h, R = -9(-y^)9(y^ - y^J^^r (1 - e -W* ) , (A-12) 

Note that the central-cut diagram in Fig. 6 corresponds to the interference between t and 
w-channel of the qq — > gg annihilation processes in Fig. 5. Since the splitting function is 
symmetric in z and 1 — z, a factor of 2 comes from the fragmentation of both gluons in 
the central-cut diagram. 

The s-channel of qq — > gg is shown in Fig. 7 which has only one central-cut. Its contri- 
bution to the partonic hard part is, 



—d a 2 s x B [df T }dz- _ , , N 2(z 2 - z + l) 2 C F C A 



H "> = if J if J t' 7 - c \ {1 - z) - -jjf ' (A " 13) 



I 7 C = 0(-y-)6(y~ - y^)e' (x+XL)p y e~ lXLP {y ^ 'e _t!e£P . (A-14) 

Note that the splitting function 2(z 2 - z + l) 2 /z(l - z) = 2[1 - z(l - z)] 2 /z(l - z) is 
symmetric in z and 1 — z. Therefore, fragmentation of the two final gluons gives rise to 
the factor of 2 in front of the gluon fragmentation function. 



34 



Fig. 8. The interference between t and s-channel of qq 



gg annihilation. 




Fig. 9. The complex conjugate of Fig. 8. 

The interferences between t and s-channel of qq —> gg processes are shown in Figs. 8 and 
9. There are only two possible cuts in these diagrams. The contributions from Fig. 8 are: 
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(A-15) 



(A-16) 

(A-17) 
(A-18) 
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Contributions from Fig. 9, which are just the complex conjugate of Fig. 8, are: 
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CfQa 
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(A-19) 



(A-20) 

(A-21) 
(A-22) 



One can collect all contributions of the double hard qq — > gg processes from the central- 
cut diagrams, which should have the common phase factor 



i c = e(- y -)6{y~ - y-) e ix p + y~ e - iXLP+{y z } 



(A-23) 



and obtain the total effective splitting function in the hard partonic part, 



C 



2 1 



N P q ^gg{z) z{l _ z)Nc 



{C 2 F [1 + z 2 + 1 + (1 - z) 2 } - 2C F {C F - C A /2) 



+ 2C F C A {1 -z + z 2 ) 2 - C F C A [l + z 3 + 1 + (1 - z) 3 }} 
C F 



2C / 2 + (1 _ " )2 - 2C A [z 2 + (1 - z) 2 ] 



(A-24) 



We will find later in Appendix A-3 that the above result can also be obtained from the 
total matrix elements squared for qq — > gg annihilation. 

We now consider the annihilation processes qq — » q^ with qi ^ q. There is only the 
s-channel process with one central-cut diagram as shown in Fig. 10. Its contribution to 
the hard part is 
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Fig. 10. s-channel qq — > qiqi annihilation. 




Fig. 11. i-channel qqi(qi) — ► qqi(qi) scattering. 



x[z 2 + (l-zY] 



2 ,C F 



io,c 
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(A-25) 
(A-26) 



Here we define the effective splitting function for qq — ■> annihilation as, 
^W(*) = §V + (1-*) 2 ]- 



iV, 



(A-27) 



Similarly, for ggj — > qqi scattering with qi ^ q, there is only the t-channel as shown in 
Fig. 11. There are, however, three cut diagrams. Their contributions to the partonic hard 
part are: 



—d a 2 s x B f dl\ f dz- 
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, , v l + Z 2 , . 1 + (l-z) 2 
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(A-2* 
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(A-29) 

(A-30) 
(A-31) 
(A-32) 



The twist-four two-parton correlation matrix element associated with the above quark- 
antiquark scattering is the quark- ant iquark correlator, 



T4 (x, x l ) oc e ^P + y-- i ^P + (v2 -vi) 

x(A|^(0)^(y-)^(yr) , (A-33) 

and one should sum over all possible q%^ q flavors. Note that in the above matrix element, 
the momentum flow for the antiquark (qi) is opposite to that of the quark (q) fields. 

For quark-quark scattering, qqi — > qq^ the hard part is essentially the same. The only 
difference is the associated matrix element for the quark-quark correlator which is ob- 
tained from that of the quark-antiquark correlator via the exchange i\) qi {y<i) —>■ ^(2/2) 
and tpqXyi) -> Mvi), 



T*.(x, x l ) oc e ixp+y +ix ^ + ^ ) 



X{A\^M\^{V~)MV2) \4>M\A) ■ (A-34) 

Note that the momentum flows of the two quarks (q and qi) point in the same direction. 

The effective splitting function of this scattering process is defined through the fragmen- 
tation of the quark in the central-cut diagram, 

Cf p / \ Cp 1 + z 2 



For annihilation qq — > qq into identical quark and antiquark pairs, in addition to the 
s-channel (Fig. 10 for qi = q) and t-channel (Fig. 11 for q^ = q), one has also to consider 
the interference between s and t-channel amplitudes as shown in Figs. 12 and 13, each 
having two cuts. Their contributions to the hard partonic parts are, respectively: 
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Fig. 12. Interference between s and ^-channel of qq — > qq scattering 




Fig. 13. The complex conjugate of Fig. 12. 
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(A-37) 

(A-38) 
(A-39) 
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Fig. 14. The interference between t and it-channel of identical quark-quark scattering qq — > qq. 
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One can again collect contributions from the central-cut diagrams of the double scattering 
processes in Figs. 10, 11 12 and 13 and obtain the total effective splitting function for 
qq -> qq, 



^P- -I 
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z z + (l 



C F 1 + z 2 _ C F (C F -C A /2) Az< 
l + z 2 1 2z 2 

T. ^ + 



N r 
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Here we have used C F — Ca/2 = —1/2N C . For antiquark fragmentation, P c 
Pqq-+qq(l — z) . One can also obtain the above result from qq 
squared as shown in Appendix A- 3. 



qq^qq\ 



qq scattering matrix 



Similarly, for scattering of identical quarks qq — > qq, one should set q^ = q in Fig. 11 [in 
Eq. (A-28)]. In addition, one should also also include interference between t and w-channel 
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of the scattering as shown in Fig. 14. The contributions from such interference diagram 
are, 



-jj-D _ a 2 s x B f dt\ r dz 7 
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(A-46) 

(A-47) 
(A-48) 
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Note again that the fragmentation of both quarks contributes to the factor 2 in Eq. (A- 
45) since the splitting function is symmetric in z and 1 — z. The twist-four two-quark 
correlation matrix element associated with qq —> qq scattering is compared 
to T^(x,xl) for quark-antiquark annihilation processes. 

We can sum contributions from the double hard scattering in all the central-cut diagrams 
in Figs. 11 and 14 and obtain the total effective splitting function for qq — > qq processes, 
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There are two remaining cut diagrams that contribute to the quark-antiquark annihilation 
at the order of 0(a 2 ) as shown in Figs. 15 and 16. Their contributions are: 
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Fig. 15. Interference between final-state gluon radiation from single and triple-quark scattering. 
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Fig. 16. The complex conjugate of Fig. 15. 
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A-2 Effective splitting functions 



In this Appendix, we list the effective splitting functions associated with each process 
qa — > b and the double-hard {HI), hard-soft (SI) or their interferences (1,12) according 
to Eq. (73). 
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The non-singlet splitting functions for qq — > b, defined as 
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are listed as below: 



AP. 



(/) (z) 
qq^qi(qi) { >' 



qq^qM)( Z "> ~ P qq^q*(q*)( Z ^ ^ P qq^qr{qi)^ ~ P <, 

(1 - Z 2 ){1 + Z 2 + (1 - z) 2 ) 2 1 + Z* 



Ap N(HI) 



» 



qq^q 

K z) = P c 



jN{HI). 



N c z(l-z)' 



, z) AP 
qq^qy^li L -^ ± qq-^g 



z)=P c 



{HI) (z) 



43 



2C, 



l + z' 1 + (1 



4«(«) 



2C F i±^ + 2 



l + (l-z) 



1-Z N r z(l-Z) 



P^Uz), AP^(z) = P^) (6 = g,g,^) 



,JV(/2) , 
qq- 



(A-65) 

(A-66) 



A-3 Alternative calculations of central-cut diagrams 



As a cross-check of the hard partonic parts calculated from different cut diagrams in 
Appendix A-l, we provide an alternative calculation of all the central-cut diagrams, 
which correspond to quark-quark (antiquark) scattering. 

Considering a parton (a) with momentum q scattering with another parton (b) that 
carries a fractional momentum xp, a(q) + b(xp) — > c(£) + d(p'), the cross section can be 
written as 



o 4 d 3 £ 
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where g = [0, g ,0] and p = 0, 0] are momenta of the initial partons and 



f2 
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2zq~ 



zq~ 
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is the momentum of one of the final partons. With the given kinematics, the on-shell 
condition in the cross section can be recast as 



(xp + q — £) 2 = 2(1 — z)xp + q [l ) 

The Mandelstam variables of the collision are 



2z(l — z)p + q~ 
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where we have used the on-shell condition x = xl- 

With Eq. (A-67) and parton distribution functions f b N (x), one can obtain the parton- 
nucleon cross section, 
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where s = 2p + q is the center-of-mass energy for aN collision, Co is some common color 
factor in the scattering matrix elements and 



P^cd(z) = (l/C )\M\ 2 ab _> cd (t/s,u/s) 
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is what we have defined as the effective splitting function for the corresponding processes. 
One can therefore easily obtain these effective splitting functions from the corresponding 
matrix elements for elementary parton-parton scattering [39]. We will list them in the fol- 
lowing. A common color factor for all quark-quark(antiquark) scattering is Co = Cf/N c . 

qq — > qiCji annihilation: 
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qq — > qq annihilation: 
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qq — > gg annihilation: 
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?9t(9i) scattering: 
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qq — > gg scattering: 
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For quark-gluon Compton scattering, the relevant gluon distribution function is xlGn(xl) 
One can therefore rewrite contribution from qg — > qg to Eq. (A-71) as, 
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We have then for qg — ► qg scattering, 
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Comparing this result with that in Ref. [18] for the quark-gluon rescattering, we can see 
that they agree in the limit 1 — z — > 0. This is a consequence of the collinear approxi- 
mation employed in Ref. [18] in the calculation of the hard partonic part in quark-gluon 
rescattering. 



We can also extend this calculation to the case of gluon-nucleon scattering. One can use 
Eq. (A-71) to define the splitting function for gq — > gq scattering, 
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Here for gluon-parton scattering, there is no common color factor. 

99 ~ > <7<7 annihilation, 
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99 ~^ 99 scattering 
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One can use this technique to extend the study of modified fragmentation functions to 
propagating gluons. Since the modification is dominated by quark-gluon and gluon-gluon 
scattering, comparing the effective splitting functions, 
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in the limit z — > 1, one can conclude that a gluon's radiative energy loss is larger than 
a quark by a factor of N C /C F = Ca/C f = 9/4. We will leave the complete derivation of 
medium modification of gluon fragmentations to a future publication. 
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